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ON THE FOURIER ANALYTIC STRUCTURE 
OF THE BROWNIAN GRAPH 


JONATHAN M. FRASER AND TUOMAS SAHLSTEN 

Abstract. In a previous article {Int. Math. Res. Not. 2014, 2730-2745) T. Orponen and 
the authors proved that the Fourier dimension of the graph of any real-valued function on M 
is bounded above by 1. This partially answered a question of Kahane (’93) by showing that 
the graph of the Wiener process Wt (Brownian motion) is almost surely not a Salem set. In 
this article we complement this result by showing that the Fourier dimension of the graph 
of Wt is almost surely 1. In the proof we introduce a new method based on ltd calculus to 
estimate Fourier transforms by reformulating the question in the language of ltd drift-diffusion 
processes and combine it with the classical work of Kahane on Brownian images. We also give 
some applications to the equidistribution of orbits under toral endomorphisms. 


1. Introduction and results 

1.1. Geometric properties of Brownian motion. Gaussian processes are standard models 
in modern probability theory and perhaps the most well studied example is the Wiener process 
(or standard Brownian motion) W = Wt dehned such that Wq = 0, the map t ^ Wt is almost 
surely continuous and Wt has independent increments such that Wt — Wg for t > s is normally 
distributed; 

Wt-Wg^ N{0,t-s). 

The Wiener process has far reaching importance throughout mathematics and it is a topic of 
particular interest to understand the geometric structure of the process. This can be achieved 
by studying several random fractals associated to the process such as images W(K) := {Wt : 
t G K} of compact sets K C [0,cx)), level sets Lc(W) := {t G M : hU = c} for c G M, graphs 
G{W) := {(t, Wt) : t G M} and other more delicate constructions such as SLE^-curves. 

Due to the properties of the Brownian motion, these random fractals enjoy a certain statistical 
self-similarity which facilitates computation of their Hausdorff dimensions dimp. Classical 
results include McKean’s proof [27] that dimHlU(A) = min{l, 2 dimp almost surely for 
each compact K C [0, cx)). Moreover, for the level sets dimHTc(lU) = 1/2 almost surely for 
c = 0 by Taylor [37] and for all c G M by Perkins [31] whenever L^iW) is non-degenerate. 
For the Brownian graph G{W), Taylor [36] proved that dimHG(lU) = 3/2 almost surely and 
Beffara computed Hausdorff dimensions of SLE^-curves [3]. Moreover, Hausdorff dimensions 
for similar sets given by many other Gaussian processes, such as fractional Brownian motion. 
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have been also considered, see for example Adler’s classical results [1] for fractional Brownian 
graphs and the recent work for variable drift by Peres and Sousi [30]. 

1.2. Fourier analytic properties of Brownian motion. The Hausdorff dimension is the 
most commonly used tool for measuring the size of a set A but there is also another fundamental 
notion based on Fourier analysis which reveals more arithmetic and geometric features of A, 
including curvature, which are not seen by the Hausdorff dimension. This is based on studying 
the Fourier coefficients of a probability measure /i on A C M'^, which are dehned by 

m--= I eeM". 

Now the size of A can be linked to the existence of probability measures fi on A with decay of 
Fourier coefficients fi{f) when |^| —)■ cx). The following connection between Hausdorff dimension 
and decay of Fourier coefficients is well-known and goes back to Salem and Kaufman, but we 
refer the reader to [25] for the details. If dimn A > s, then A supports a probability measure 
H with 1/1(01 = “on average”, that is, 1/^(0 PI'Cl* ^ df < oo and vice versa the 

Hausdorff dimension can be bounded from below if such a measure fi can be found. It is 
possible, however, that dimn A = s > 0 but no measure fi on A has Fourier decay at inhnity, 
this happens for example when A is the middle-third Cantor set on M. Therefore, one dehnes 
the notion of Fourier dimension dimp A of a set A C M'’* as the supremum of s G [0, d] for which 
there exists a probability measure fi supported on A such that 

IW{)l=0(|«|-'^), (1.1) 

Then by this dehnition we always have dimp A ^ dimn A and if the two dimensions coincide 
then A is called a Salem set or a round set after Kahane [20]. In general Fourier dimension 
and Hausdorff dimension have no relationship other than this; in fact, Korner [11] established 
that for any 0^s<t^litis possible to construct examples A C M with dimp A = s and 
dimn A = t. Further theory for Fourier dimension was recently developed by Ekstrom, Persson 
and Schmeling [8]. For a more in depth account of Fourier dimension, the reader is referred to 
[24, 25]. 

Finding measures fi on A with polynomially decaying Fourier transform (i.e. (1.1) for some 
s > 0) has deep links to absolute continuity, arithmetic and geometric structure and curvature. 
If A supports a measure fi such that (1.1) holds with s > 1, then Parseval’s identity yields that 
fi is absolutely continuous to Lebesgue measure and A must contain an interval. An application 
of Weyl’s criterion known as the Davenport-Erdos-LeVeque criterion (see [7]) yields that in M 
polynomial decay of fi guarantees that fi almost every number is normal in every base and a 
interesting result of Laba and Pramanik [13] shows that if the s in (1.1) is sufficiently close 
to 1 for a measure /i on A C M, then A contains non-trivial 3-term arithmetic progressions. 
Moreover, an analogous result also holds for higher dimensions with arithmetic patches [6]. 
On the curvature side, if A is a line-segment in then A cannot contain any measure with 
Fourier decay at inhnity so A cannot be a Salem set. However, if A is an arc of a circle or 
more generally a 1-dimensional smooth manifold with non-vanishing curvature when the 1- 
dimensional Hausdorff measure fi on A satishes (1.1) with s = 1, see [25]. In particular, A is a 
Salem set. In these examples of A one can observe that the important arithmetic or curvature 
features present are not seen from the Hausdorff dimension. 
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Constructing explicit Salem sets (which are not manifolds), or just sets A supporting a 
measure /r satisfying (1.1) for some s > 0, can be achieved through, for example, Diophantine 
approximation by Kaufman’s works [22, 23], Bluhm [4] or via thermodynamical tools by Jordan 
and Sahlsten [12], However, for random sets it has been observed in many instances that A is 
either almost surely Salem or at least supports a measure /i with (1.1) for some s > 0. This was 
first done for random Cantor sets by Salem [34], where Salem sets were also introduced. Later 
Kahane published his classical papers [16, 17], where he found out that the Wiener process and 
other Gaussian processes provide natural examples. 

The study of Fourier analytic properties of natural sets for Gaussian processes and more 
general random fields has since been an active topic. For the Brownian images Kahane [19] 
proved that for any compact iC C M the image W{K) is almost surely a Salem set of Hausdorff 
dimension min{l, 2 dimiC}. Kahane also established a similar result for fractional Brownian 
motion. Laba and Pramanik [13] then applied these to the additive structure of Brownian 
images. Later Shieh and Xiao [35] extended Kahane’s work to very general classes of Gaussian 
random fields. However, understanding the Fourier analytic properties of the level sets and 
graphs remained an important problem for some time. In 1993, Kahane [20] outlined the 
problem explicitly. 

Problem 1.1 (Kahane). Are the graph and level sets of a stochastic process such as fractional 
Brownian motion Salem sets? 

In this form Kahane’s problem was formulated by Shieh and Xiao [35, Question 2.15]. For the 
Wiener process Kahane [18] had already established that the level sets Lc{W) are Salem almost 
surely for any fixed c G M provided that Lc{W) is non-degenerate. The fractional Brownian 
motion case has recently been considered in the case c = 0 by Fouche and Mukeru [9]. 



Figure 1. Three realisations of the graph G{W) for the Brownian motion Wt- 


Kahane’s problem for graphs, even in the case of the standard Brownian motion IW, however, 
remained open for quite a while until, together with T. Orponen, we established that the 
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Brownian graph G{W) is almost surely not a Salem set [14], It turned out that the reason for 
this is purely geometric: the proof was based on the following application of a Fourier-analytic 
version of Marstrand’s slicing lemma. 

Theorem 1.2 (Theorem 1.2 in [14]). For any function f : [0,1] —)■ M the Fourier dimension of 
the graph G{f) cannot exceed 1. 

Indeed, since the Hausdorff dimension dimHG(IF) = 3/2 > 1 almost surely (see [36]), this 
result shows that the graph G(W) cannot be a Salem set almost surely, which answers Problem 
1.1 in the negative for the Wiener process. Note that this also gives a negative answer for 
fractional Brownian motion since the Hausdorff dimension in that case is also strictly larger 
than 1 almost surely. 

The methods in [14] are purely geometric and involve no stochastics, so they do not shed 
any light on the precise value for the Fourier dimension of G{W). Note that even though 
dimHG(/) ^ 1 for any continuous / : [0,1] —)■ M by Theorem 1.2, it is completely possible that 
the Fourier dimension dimFG(/) = 0. For example this happens by just taking / to be affine 
or even any Baire generic f G C[0,1] satishes dimp G(/) = 0, see [14, Theorem 1.3]. 

The main result of this paper is to complete the work initiated by Kahane’s Problem 1.1 
in the case of Brownian motion by establishing the precise almost sure value of the Fourier 
dimension of G{W). 

Theorem 1.3. The graph G{W) has Fourier dimension 1 almost surely. 

Moreover, the random measure p we use to realise the Fourier dimension is Lebesgue measure 
dt on [0,1] lifted onto the graph G{W) via the mapping t h-)■ (t, Wt). The precise estimate we 
obtain is that almost surely 

IKOl = 0(|^rVlog|^l)> asl^l^oo, (1.2) 

which combined with Theorem 1.2 yields Theorem 1.3. 

A natural direction in which to continue this line of research would be to study other Gaussian 
processes with different covariance structure, see the discussion in Section 5.1. 


1.3. Methods: Ito calculus and reduction to Brownian images. The key method we 
introduce to estimate the Fourier transform of the graph measure p, is based on ltd calculus, 
which has previously been a natural framework in the theory of stochastic differential equations. 
As far as we know, Ito calculus has not been previously considered in this Fourier analytic 
context. In Section 5 we will discuss the prospects of this approach. Next we will discuss 
this method and give a brief summary of the main steps in the proof. When written in polar 
coordinates, (1.2) asks about the rate of decay for the integral 

[ exp{—27riu{tcos9 + Wtsm9))dt 
0 

for f = u{cos9,sm9) G n > 0, 6* G [0,27r), as m —)■ cx). There are two distinct cases we will 
consider depending on the direction of f, which we give a heuristic description of here. 

If we ignore the random component Wt sin 9, that is, set 0 = 0 or tt, then standard integration 
using the chain rule shows that fi{^) equals the Fourier transform of Lebesgue measure dt at 
u, which decays to 0 with the polynomial rate u~^ = so we are done for these directions. 
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However, if 6 is not equal to 0 or tt, we still have a small random (non-smooth) term hhtsin^, 
so a classical change of variable formula or other tools from classical analysis cannot be used. 

The new approach we introduce is to utilise tools from stochastic analysis. The key obser¬ 
vation is that we can write = f exp(iXt) dt, where the stochastic process 

Xt := bt + aWt 

is a so called ltd drift-dijfusion process, where b := —27rMcos6* is the drift coefficient of Xf 
and a := —27rasin0 is the diffusion coefficient of Xf. ltd processes have many useful analytic 
tools from ltd calculus (see Section 2) associated to them, in particular ltd’s lemma, which 
works as an analogue for the chain rule. The price we pay is that ltd’s lemma introduces some 
multiplicative error terms involving stochastic integrals, but they can be estimated with other 
tools from ltd calculus using moment analysis. 

The estimates we obtain from ltd calculus allow us to obtain the correct Fourier decay (1.2) 
for /i when 9 is close to 0 or tt with respect to u~^ (more precisely, |sin0| < In 

other words, when f is close to pointing in the horizontal directions. Thus another estimate 
is needed for 6 bounded away from 0 and tt. This is where Kahane’s classical work [19] on 
Brownian images comes into play. If we completely ignore the deterministic component tcos6, 
by setting 6 = 7 r /2 or 37 r/ 2 , then fi{f) is the Fourier transform the Brownian image measure 
u, that is the t ^ Wt push-forward of the Lebesgue measure dt on [0,1] at u. Kahane [19] 
in fact already established that almost surely the decay of \i'{u)\ is almost surely of the order 
u~^^/\ogu = |,^|“^A/log |.^| so (1.2) holds for these directions. Now a simple modihcation to 
Kahane’s argument reveals that whenever 6* 7 ^ 0 or tt, then almost surely 

1/1(01 = 0(1 sin0|-^|O“Vlog|^l) 

see the discussion in Section 3.3. Now one notices that when 6 approaches 0 or tt, this estimate 
blows up, and so one cannot obtain a uniform estimate over all directions from this. However, 
this gives ( 1 . 2 ) if |sin 0 | ^ so combining with the estimates we obtained through ltd 

calculus, we are done. See Section 3 for more details on the main steps of the proof. 

1.4. Equidistribution. As mentioned in the introduction, Fourier coefficients of /i on a set A 
are deeply linked to arithmetic properties of the set A. In this section we will point out one 
application one can achieve from Theorem 1.3 (in particular the decay (1.2)). Let T := M/Z 
be the unit circle or 1-torus, which is naturally identihed with [0,1]. A sequence (xn) C M is 
equidistributed if 

tj{l ^ n ^ A^ : Tn e [a,b]] 

N 

as A^ —)■ cx) for any interval [a, b] C [0,1]. If the sequence Xn = p^'x mod 1 equidistributes for 
a hxed p E N, then we say x is normal in base p. Checking normality of a specihc x is often 
difficult so instead one often tries to check if a random number (with respect to a measure p) 
is normal. This is where the decay of the Fourier coefficients of p are crucial: 

Theorem 1.4 (Davenport-Erdos-LeVeque). Let p be a measure on M for which there is some 
a > 0 such that 


Ia^(0I = 0(1^1 “), as 1^1 ^ 00 . 
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Then for any strictly increasing sequence of reals (<s„)neN the sequence {snX mod 1) is equidis- 
tributed for fi almost every x. In particular, p almost every x is normal in every base. 


Davenport-Erdos-LeVeque criterion was established in [7] and in this form by Queffelec- 
Ramare [33] and it is a consequence of the classical Weyl’s equidistribution criterion (1908), 
which says that a sequence (xn) C T equidistributes if and only if the exponential sums 

N 

I 'y^^exp{27rikxn) = o{N), as iV —)■ cx), 


71=1 


for any k &'L \ {0}. In higher dimensions one can also study analogues of equidistribution on 
the d-torus := /!/'. Here a sequence {xn) C is equidistributed if in the weak topology 

1 ^ 


lim — = Haar measure on T'^. 

JV^OO N ^ " 


n=l 


Weyl’s criterion in the higher dimensions has the following form; 

Theorem 1.5 (Weyl’s equidistribution criterion). A sequence (xn) C is equidistributed if 
and only if for any k G \ {0} we have 


1 


N 


lim — exp(27rf k ■ Xn) = 0. 
N^oo N ^ ^ ’ 


71=1 


The natural analogue of sequences [xn) on the torus is to consider orbits under suitable 
expanding maps. For example, being normal is base p precisely means that for the xp map 

Tp{x) := px mod 1 

the orbit {T”a;}neN of x is equidistributed. The map Tp is an example of an expanding toral 
endomorphism, that is, p'^x grows in norm as n —)■ cx and Tp preserves the algebraic structure 
of the unit circle T. A similar notion can be dehned for toral endomorphisms on T'’*, that is, 
maps of the form 

T{x) = Ax mod 1 

for some matrix A G GLrf(Z). We say that the map T is expanding if the singular values of 
A are all strictly greater than 1. We prove the following higher dimensional analogue of the 
Davenport-Erdos-LeVeque criterion; 

Theorem 1.6. Let p be a measure on for which there is some a > 0 such that 

l/^(OI = 0(|^r“), as 1^1 ^x. 

Then the T-orbit {T"(a;)}„gN of x equidistributes at p, almost every x for any expanding toral 
endomorphism T ; —)■ T^. 


The proof of Theorem 1.6 follows the same lines as the M-version by Davenport-Erdos- 
LeVeque and we postpone the proof until Section 4. Thus as an immediate corollary of Theorem 
1.3, we obtain the following equidistribution theorem for Brownian motion; 

Corollary 1.7. Almost surely for Lebesque almost every t G [0,1] the orbit {T'^{t,Wt)}n£f>i 
equidistributes for any expanding toral endomorphism T ; —)■ . 
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Figure 2. A realisation of the graph of Brownian motion Wt embedded onto 
the 2-torus with the beginning of a T-orbit oi x = {t, Wt). 

1.5. Organisation of the paper. In Section 2 we give the necessary background from ltd 
calculus. In Section 3 we will give the proof of our main result Theorem 1.3. The key estimates 
are obtained in Section 3.2 and Section 3.3 corresponding to the two cases discussed above. In 
Section 4 we prove the equidistribution Theorem 1.6 and finally in Section 5 we discuss further 
prospects for this work. 


2. ItO CALCULUS 

2.1. Stochastic integration. In the proof of the main Theorem 1.3 we will end up studying 
integrals of the form J f{Xt) dt for some stochastic processes Xt and smooth scalar functions 
/. As standard analysis methods cannot be applied to these integrals, we need theory from 
stochastic analysis. Stochastic analysis provides a pleasant framework to deal with non-smooth 
processes, like the Wiener process Wt, but still preserves much of the analytic tools available 
in smooth analysis. In this section we discuss the specihc tools from ltd calculus which we will 
rely on, the main references for this section are given in the book [21]. 

Let (n, X, {Xt)t^o, P) be a filtered probability space, that is, Xt C X is an increasing filtration 
in t. Let VF = W be the Wiener process adapted to this filtered probability space, that is, 
Wt is Xt measurable and for each t,s ^ 0 the increment Wt+s ~ Wt is independent of Xt- We 
say that an M or C valued stochastic process Zt is adapted if it is Xt measurable for all f ^ 0. 
We will say that a real or complex valued adapted process Zt is Wt integrable if the quadratic 
variation \Zt\^ dt < oo for any time T ^ 0. Given a real valued adapted Wt integrable 
stochastic process Xt, then P almost surely for any time T ^ 0 it is possible to construct a 
stochastic integral 

XtdWt 

of Xt with respect to Wt in the sense of ltd, see [21, Chapter 3.2]. We use the differential 
notation dUt = Xt dWt to mean that P almost surely Ut — Uq is the stochastic integral of Xt 
with respect to Wt at time T ^ 0. 
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We mainly deal with complex valued stochastic processes so for the sake of convenience we 
will also dehne the complex valued stochastic integral for a C valued Wt integrable adapted 
process Zt is dehned coordinate-wise using real integrals: 

rT rT rT 


ZtdWt:= / ReZtdWt + t 


Im Zt dWt, 


Jo Jo Jo 

where the real integrals are standard M valued stochastic integrals with respect to the Wiener 
process Wt. We write dZt := dXt + idYt for a complex valued process Zt = Xt + iYt with M 
valued Xt and Yt. 


2.2. Ito drift-diffusion processes. The main class of adapted processes where ltd calculus 
is applied to and which will appear in our proof of Theorem 1.3 are given by Wiener processes 
with drift and diffusion coefficients. These are called Ito drift-diffusions; 


Definition 2.1 (Ito drift-diffusion process). A real or complex valued adapted stochastic pro¬ 
cess Xt is called a ltd drift-dijfusion process if there exists a Lebesgue integrable adapted bt 
and Wt integrable adapted at such that Xt satishes the stochastic differential equation 

dXt = btdt + at dWt- 

For Ito drift-diffusion processes there exists the following important analogue of the change of 
variable formula, which follows from robustness of Taylor expansions for stochastic differentials: 

Lemma 2.2 (Ito’s lemma). Let Xt he an ltd drift-dijfusion process and / : M —)■ M twice 
differentiable. Then /(W) is an ltd drift-dijfusion process such that P almost surely for any 
T ^ 0 we have 

fiXr) - /(Xo) = £ (btfiXt) + ^/"(W)) dt + £^tf{Xt) dWt. 

Ito’s lemma in this form was given in this pathwise form in [21, Theorem 3.3]. By using 
the dehnition of complex valued stochastic integral, we can also write a complex valued Ito’s 
lemma as follows: 


Lemma 2.3 (Complex Ito’s lemma). Let Xt be a ltd drift-diffusion process and / : M —)■ C 
twice differentiable. Then f{Xt) is an ltd drift-dijfusion process such that for P almost surely 
for any T ^ 0 we have 

fiXr) - /(Xo) = £ (btffXt) + ^nXt)) dt + £<^tf{Xt) dWt. 

Proof. We can write / = /i + if 2 for real valued twice differentiable /i ,/2 : M —)■ M. Then the 
derivatives f' = f[+ i /2 and f" = /" -1- iff. Moreover, by Ito’s lemma (Lemma 2.2) we obtain 
for each j = 1,2 that 

2 

df,{Xt) = (btf'iXt) + ^f'iXt)) dt + atf'iXt) dWt. 

Then by the convention df{Xt) = dfi{Xt) + idf 2 {Xt) this gives 

2 

df{Xt) = (btfiXt) + ^nxt)) dt + atfiXt) dWt 
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as required. □ 

2.3. Moment estimation. Ito’s lemma allows us to pass from integrals of the form f{Xt) dt 
to fj' g(Xi) dWt for functions g obtained from derivatives of /. In our case we will end up trying 

to understand the higher moments of the stochastic integrals g(Xt) dWt, which will tell us 
about the distribution of these integrals. A very standard tool to compute the moments in ltd 
calculus are the ltd isometry and more general Burkholder-Davis-Gundy inequalities (see [5]), 
which allows to pass from stochastic integrals to their quadratic variations (that just involve 
Lebesgue integral). Here we state the following version; 


Lemma 2.4 (Burkholder-Davis-Gundy inequality). Let Xt be a real valued Wt integrable 
adapted process. Then for all 1 ^ p < oo vue have 


E 


sup 


X,dW, 


2p-i 


^ 2a/10p E 



This was given by Peskir [32] with the constant 2^10p. 


3. Proof of the main result 


3.1. Preliminaries and overview of the proof. Let us now review how we will prove (1.2) 
and thus Theorem 1.3. Fix ^ = u{cos6,smd) G with modulus u > 0 and argument 
6 G [0,27r). Notice that by the definition of the graph measure p, the Fourier transform has 
the form ^ 

/i(0 = / exp(fXi) dt, 

Jo 

where Xt is the real valued stochastic process 

Xt '.= —27ru{t cos 9 Wt sin 9). (3.1) 


The first observation is that Xt is an adapted Wt integrable process and in fact an Ito drift- 
diffusion process (recall Definition 2.1) satisfying 

dXt = bdt -\- a dWt 


for deterministic and time independent coefficients b = —2tiucos9 and a = —2Tiusin9. The 
proof of bounding fl{^) will heavily depend on the value of the angle 9 we have for ^ and in 
particular how close the determining angle 9 is to 0, tt or 2n is with respect to For this 

purpose, we define the notions of horizontal and vertical angles: 


Definition 3.1 (Horizontal and vertical angles). Define the threshold angle 

9u := min{u"^/^ f}. 

Partition the angles [0,27r) using 9u into the horizontal angles 

Hu := [0, 9u] U[7r - 9u,7r + 9u] U [27r - 9u, 27r). 

and the vertical angles 

Vu := [0,27r) \ Hu- 

In other words Hu contains the 9u neighborhoods of 0 and tt on the circle mod 27r and 14 fhe 
7r/2 — 9u neighborhoods of n/2 and 37r/2 respectively, see Figure 3.1. 
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TT 



Stt 

2 


Figure 3. Splitting of [0, 27r) to horizontal angles and the vertical angles 14 

The proof will split into two cases in Sections 3.2 and 3.3 for bounding the Fourier transform 
depending on whether 6 G or 6^ G 14 - 

(1) Section 3.2 concerns angles 9 G FTn, that is, close to horizontal directions 0 or tt and as 
mentioned in the introduction our main hope here is that the smallness (with respect to 

of the diffusion component hWt will help us in transferring the decay of Lebesgue 
measure to the decay of Ji. This is where in particular Ito’s lemma (see Lemma 2.3) 
becomes crucial as it can be applied to the process f{Xt) with the function f{x) = 
exp(ix). 

(2) Section 3.3 handles the angles 6 & Vu and here the plan is to use the fact that we 
are bounded away from horizontal angles to ignore the drift component bt of the 
drift-diffusion process Xt and apply Kahane’s bound for these directions. This turns 
out to be possible due to a representation of the higher moments Kahane obtained in 
his result on Brownian images. 

It turns out that in both of the Sections 3.2 and 3.3 we only obtain decay of the Fourier 
transform /i(k) for k in an e-grid el? for all small e > 0. Here the randomness will depend 
on e > 0 but thanks to an argument also used by Kahane in [19], one can pass from this 
information to the full decay almost surely. See Section 3.4 for the details. 

Let us now proceed to bound |4(0I- both Sections 3.2 and 3.3 below we will end up 
bounding trigonometric functions with respect to 6u and for this purpose we will need the 
following standard bounds: 

Lemma 3.2 (Trigonometric bounds). ITe have the following bounds: 

(1) If 9 ^ Hu, then 

I sin 6^1 4 and jcos^l^:^- 
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(2) Ifee Vu, then 

I sin6'| ^ 

Proof. For a G [0, |] we have that both cos and sin are non-negative. Moreover, here ^ 
sin a ^ a. Thus for 9 G [0,6'„] we have 

sin6* ^ 0 ^ ^ and cos6' ^ cos6'„ ^ cos | 

and for 9 G {9u, f ] as sin | we obtain 

sin6' ^ 

This gives the claim as we may reduce the estimates back to the estimates for 9 G [0, |] by 
using standard invariance identities for sin and cos. □ 


3.2. Horizontal angles. When 9 G Hu we will first obtain the following estimate on e-grids: 


Lemma 3.3. Fix e > 0. Almost surely there exists a random constant Cuj > 0 such that for 
any k = u{cos9,sm9) G \ {0} with 9 G Hu we have 


Given ^ = u{cos9, sin0) G \ {0} and a realisation (Wt) dehne a random time T = ^ 

[0,1], which the minimum value of t G [0,1] such that 




—27r[n(cos6'-|-IFi sin6*)], if Xi ^ 0; 
—27r[M(cos6'-|-IFi sin6*)J, if Xi < 0. 


Such a time T exists almost surely since Xq = 0 and Xt is almost surely continuous (since Wt 
is almost surely continuous). Splitting the integral of Zt up into ‘complete rotations’ and ‘what 
is left over’, one obtains 


Zt dt = 


Ztdt + / Zt dt. 


For the integral over [T, 1] we get the following estimate. 


Lemma 3.4. Almost surely there exists a random constant > 0 such that for any ^ = 
u{cos9,sm9) G \ {0} with 9 G Hu we have 

, ri 


Zf dt 




- 1/2 


Proof. Since Wt is almost surely continuous, there almost surely exists a random constant 
Mu, > 1 such that Wt G [—Mu,, MJ\ for all t G [0,1]. Dehne the real-valued process 

Yt := u{t cos 9 + Wt sin 9) 

so Xt = —27iYt. Suppose Xi ^ 0. In this case Yt = [Fi] ^ 0 and so W -|- 1 ^ Ft ^ Fi. 
Moreover, when Xi < 0 we have Ft = [FiJ > 0 and Fi ^ Ft ^ Fi — 1. Thus no matter what 
is the sign of Xi is, we always have almost surely 


m(cos 6* -|- hFi sin 6^) -|- 1 ^ u(T cos 9 + Wt sin 9) ^ m(cos 6^ -|- hFi sin 9) — 1. 
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Therefore, in the case cos 6 > 0 we obtain 


^ , ,,, sin6* sin6' 1 

T ^ 1 + Wi - - - Wt -^-^ 

cos U cos U u cos U 


and when cos 6^ < 0 we have 


^ , ,,, sin6' sin6* 1 

T ^ 1 + hhi- - — Wt -^ H- 

cos U cos U u cos U 


Since u G Hu Lemma 3.2 together with Wt G [— MJ\ yields 

72 


T ^ 1 - 2V2M^u-^^^ - 

u 


Recalling M^, > 1 this gives 


Zt dt 


^ / \Zt\dt = l- T 


as required. □ 

We now estimate the integral over [0,T], which is where ltd calculus comes into play. 

Lemma 3.5. Fix £ > 0. Almost surely there exists a random eonstant Cai > 0 sueh that for 
any k = u{cos9,sm9) G \ {0} with 9 G Hu we have 

rT 


Zt dt 




To prove Lemma 3.5, we hrst need to compute the higher order moments of the random 
variable Zt dt. 

Lemma 3.6. For any p G N and ^ = u{cos9, sin6') G \ {0} with 9 G Hu the 2pth moment 

|2p 


E 


Zt dt 




Proof. Recall that 

Xt = —27iu{t cos 9 + Wt sin 9) 

is an ltd drift-diffusion process satisfying the stochastic differential equation 

dXt = bdt + a dWt 


for deterministic and time independent coefficients b = —2'kucos9 and a = —2'Kusm.9. Writing 

f{x) := exp(ia;), a; G M, 

then Zt = f{Xt), f'{x) = ieyrpiix) and f'\x) = — exp(ia;). Thus by complex ltd’s lemma (see 
Lemma 2.3) we have P almost surely 

^ f{Xt)dt + ai r f{Xt)dWt. (3.2) 

Jo 

Note that ^ 1 is random and only J^i measurable, thus it is not a stopping time. However, 
as Lemma 2.3 is given pathwise, that is, P almost surely ltd’s lemma holds for any time T ^ 0 


f{XT) - /(Xo) = {U - 7/2) [ 

Jo 
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then as is P almost surely well-defined, we have (3.2) almost surely. Since Xq and Xx are 
27r multiples of integers by definition, we have /{Xx) = /(df^o) = 1- Thus (3.2) gives 


f{Xt)dt = 


ai 


hi - Jo 


f{Xt)dWt 


Since b and a are deterministic, this yields that the 2pth moment 


E 


f{Xt)dt 


2p 


ai 


2p 


E 


f{Xt)dWt 


2p 


bi — I2 I 

Applying the Burkholder-Davis-Gundy inequality (see Lemma 2.4) for the process cosX^ gives 

2p-| 


E 


r r 

2p- 


Y 

/ cosXtdH’t 


< E 

I sup 

- Jo 

- 


X 


cos Xt dWt 


^ 2 a / 10 p j cos^ Xtdt^ 

^ 2^/mp 


since cos^ ^ 1. Similar application for the process sinX* gives 


E 


sin Xf dWf 


2p- 




By Euler’s formula, we can write f{Xt) = cosXt -|- f sinX* and so 


f{Xt)dWt= / cos XtdWt + i / sin XtdWf 


Hence 


E 


f{Xt)dWt 


2p 


= E 


cos Xf dWf 


2 


2\ p-| 

+ 

/ sinX^dlTt 

) 


Jo 

/ J 



r 


2p 


2p- 

^ E 

2P 

/ cosXfdWt 
Jo 

+ 2P 

/ sinXidlTt 

Jo 

. 


= 2^* E 


cos Xf dWf 


2p- 


+ E 


sin Xf dWf 


2p-| 


^ 2^4 v^. 

Moreover, as 9 E Hu we have by Lemma 3.2 that cos^ 0^1/2 and sin^ 9 ^ u~^. Hence 


at 


a" 


sin^ 9 sin^ 9 

^ Tw = w w ^ 2u 


bi — a‘^l2\ b‘^ + a'^/A cos^ 9 cos^ 9 


Therefore, 


as required. 


E 


\2p 


f\Xt) dt ^ 44P02p ^ ^3^1/2 


□ 
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Proof of Lemma 3.5. Fix e > 0. Then for all k G eT? \ {0} define the random variable 

^(k);=(^^ Ztd1^-XAiy), 

where xa is the indicator function on the set 

A = m(cos 6', sin6') G \ {0} ; 6* G Hu}- 

Note that I(k) is well-dehned and hnite since | fj' Zt dt\ ^ 1 by | exp(ia;)| = 1 property. Lemma 
3.6 now yields for any k G \ {0} and p eN that 

E|/(k)|2p ^ 13p^/HP\k\-P 


as when k ^ A we have /(k) = 0. Write pk = L^og |k|J. Then 


E E 

keeZ 2 \{ 0 } 


1/2 ^ 
13Pk^ 4PHk|-^’>^ 


1-3 


This means that the summands tend to 0 almost surely as |k| —)■ 
random constant > 0 such that for all k G \ {0} we have 


1-3 


|/(k)| 


2pk 


13py^4Pk|k|-Pk 


^ a. 


Therefore, by possibly making bigger we obtain 


< CX). 

cx) and so we can hnd a 


i/(k)i <aikrV2, 


This holds for each k G \ {0} so by dehnition of /(k) we have whenever k 
eZ^ \ {0} with 9 G Hu that 


Zf dt 




as claimed. 


m(cos 6*, sin 6) G 


□ 


We are now in position to complete the proof of Lemma 3.3. 


Proof of Lemma 3.3. Fix e > 0. By the splitting 

Ztdt + j Zt dt 

and Lemmas 3.4 and 3.5, we have that almost surely there exists a constant > 0 such that 
for all k = M(cos0,sin0) G elZ \ {0} with 6 G Hu we have 



\m\ < 



+ 



Zt dt 


^a,lkTV2 


as required. 


□ 
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3.3. Vertical angles. In this section we apply Kahane’s work to obtain Fourier decay estimates 
when 6 & Vu- 


Lemma 3.7. Fix e > 0. Almost surely there exists a random eonstant Cai > 0 sueh that for 
any k = m(cos 6*, sin6*) G el? \ {0} with 9 E Vu we have 

|/i(k)| ^ C^|k|“^/Vlog|k|- 

Let us discuss a few estimates Kahane obtained in [19]. Let u be the push-forward of Lebesgue 
measure on [0,1] under the map t ^ Wt, that is, u is the Brownian image of Lebesgue measure. 
Kahane established the following: 

Theorem 3.8 (Kahane, page 255, [19]). Almost surely 

|z/(t;)| ^ 0(|t;|“^A/log |t|) as |t| —)■ oo. 

The key ingredient for the proof of Theorem 3.8 was based on establishing the following 
bound for the higher moments: 

Lemma 3.9 (Kahane, page 254, [19], estimate (2)). There exists a eonstant C > 0 sueh that 
for any t G M \ {0} and any p G K we have 

E|F(t;)pP ^ CPpV\v\-‘^F 

We can use Lemma 3.9 to give a bound on the higher moments in our setting, but with the 
price that the exponent will increase from —2p to —p. 


Lemma 3.10. There exists a eonstant C > 0 sueh that for any p eN and ^ = u{cos 6, sin 6) E 
\ {0} with 6 E Vu the 2pth moment satisfies 

« cnm-e 

Proof Write t = (ti,... ,tp) E [0,1]^ and dt as the Lebesgue measure on [0,1]^. Given t,s G 
[0, l]^’, we denote 

p p 

<p(t,s) := ^(4 - Sfc), 4(t,s) := ^(W^ - and T(t, s) := E|(p(t, s)p. 

k=l k=l 

By the dehnition of pe, p and the Fourier-transform, and using the fact that the multivariate 
process 

X(t, s) := —27rcos(6*)(p(t, s) — 27r sin(6*)'^(t, s) 

is Gaussian with mean —27r cos(0)(p(t, s) and variance sin^(6*)\['(t, s), we have through Fu- 
bini’s theorem and the formula for the characteristic function that 


E|/;(oU = E [ 



/ exp(—27rm(cos(6')(p(t, s) -|- sin(6*)4(t, s))) dt ds 

[o,i]p 

Eexp(mX(t, s)) dt ds 


[o,i]p 


exp(—27ri cos(6*)Mp(t, s) — 27r^|Msin(6')p\['(t,s)) dtds. 


[o,i]p 
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Thus by taking absolute values inside the integrals, and observing that |exp(ia;)| = 1 for any 
a; G M, we obtain 




exp(—27r^|n sin(6') pT(t, s)) dt ds. 


(3.3) 


a[o,i]p a[o,i]p 

On the other hand, by doing the expansion again for the Fourier transform u of the image 
measure u at v := nsin(0) G M \ {0} we see that 

E|z/(n)p^ = E / / exp(—27rin'^(t, s)) dt ds = / / exp(—27r^n^\['(t, s)) dt ds, 

J[o,i]p J[o,i]p >^[o,i]p J[o,i]p 

which equals to (3.3). Thus by Lemma 3.9 we have 

E|p(or^ ^ cpf\v\-^p. 

Since 9 G have | sind| ^ min{;^n“^/^, When | sind| ^ ^ we obtain 


F2 


□ 


CPf\v\-^P ^ {2CYfu-^P ^ {2Cffu-P. 

On the other hand, if | sind| ^ we have 

This completes the proof. 

Now we can complete the proof of Lemma 3.7 for vertical directions; 

Proof of Lemma 3.7. Fix e > 0. Then for all k = M(cosd,sind) G \ {0} dehne the random 
variable 

F(k) := /i(k)xij(k), 

where 

B-.= {i = M(cosd,sind) G \ {0} ; d G K}- 

Now F’(k) is a well-dehned hnite random variable as |/i(k)| ^ 1 for any k. From Lemma 3.10 
we obtain for any k G sT? \ {0} and p G Id that 

E|F(k)|2p ^ C^f\k\-P 

Write pk = [log |k|J. Then 

|F(k)|2Pk 


ke£Z2\{o} 


1-3 


CPkp^Pk|k|-Pk 


< 




1-3 


< oo. 


keeZ2\{0} 


This means that the summands tend to 0 almost surely as |k| —)■ cx) and so we can hnd a 
random constant Coj > 0 such that for all k G \ {0} we have 

|F(k)|2Pk 


1-3 


^ a,. 


CPkpkPk|k|-Pk 

Thus possibly making bigger this yields 

|F(k)| 
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Now this holds for each k G eT? \ {0} so by dehnition of F{k.) we have, whenever k = 
u{cos9,sm9) G \ {0} with 9 G 14, that 

|/i(k)| 4 4<^|k|“^/Vlog |k| 

as claimed. □ 

3.4. From lattices to We can now complete the proof of the main theorem. For this 
purpose, we need the following comparison lemma used by Kahane that allows to pass from 
convergence on lattices for Fourier transform to the whole space: 

Lemma 3.11 (Kahane, Lemma 1, page 252, [19]). Suppose r is a measure on with support 
in (—1,1)^. Suppose 99,4 : (0, cxd) —)■ (0, cx)) that are decreasing as t ^ 00 with the doubling 
properties 

99 (t/2) = 0{if{t)) and 'i/j{t/2) = 0{'ijj{t)) as t ^ 00 . 

If the Fourier transform of r along the integer lattice satisfies 

|r(n)| = 0((^(|n|)/4(|n|)), as |n| ^ 00 , 

then 

= C>(<^(|C|)/4(I^I)), as 1^1 ^ CX). 

Proof of Theorem 1.3. Combining Lemmas 3.7 and 3.3 we have that for any e > 0, almost 
surely, there exists some random constant > 0 such that for any k = M(cos0,sin6*) G 
£Z2 \ {0} we have 

|/i(k)| 4 Ca;|k|~^/^A/log|k|. (3.4) 

Dehne a measure on such that 

By the almost sure continuity of IW, we have that there exists a random constant iV4 > 0 such 
that the diameter of the support of p, is at most almost surely. Taking an intersection of the 
events that (3.4) holds for e = 1/n over all n G N allows us to hnd a random e = > D such 

that p is supported on a set of diameter strictly less than 1/e and (3.4) holds almost surely with 
this e. This guarantees that the measure is supported on (—1,1)^ and so applying Lemma 
3.11 with the measure r = and the maps (p{t) := y/\ogt and ^jJ{t) := t^l‘^ gives the claim. □ 


4. Proof of the equidistribution theorem 


The strategy of the proof of Theorem 1.6 is similar to that implemented by Queffelec and 
Ramare [33] for the one dimensional version of Davenport-Erdos-LeVeque. For this purpose we 
need the following Lemma. 


Lemma 4.1 (Lemme 7.3 [33]). Suppose {rjsi) C (0, cx) he a sequence of reals such that 
'12 n=i ^ Then there is a subsequence Nj —)■ x, j G K, such that 


rNj < 00 
j=i 


and 


lim 

j->oo 


Nj+1 

N, 


1 . 
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Proof of Theorem 1.6. Let /i be a measure on the d-torus for which there exists some a > 0 
and C* > 0 such that the Fourier transform satisfies 


for any ^ G \{0}. Also, let T : —)■ an expanding toral endomorphism. Fixk G Z‘^\{0} 

and for G N write 

1 ^ 

Sn{,x) := — ^exp(27rfk-T”(a;)), 

n=\ 

where k has been suppressed from the notation. By Weyl’s equidistribution criterion (Theorem 
1.5) for our claim it is enough to prove that Sn{x) —)■ 0 as —)■ cx) for p almost every x. 
Note that from this one can deduce that the result holds simultaneously for all expanding toral 
endomorphisms simultaneously follows since such maps are represented by two dimensional 
matrices with integer coefficients and so there are only countably many possibilities. 

Write 

rjy := [ \Sn{x)\‘^ dfi{x), 


and 


T{x) = Ax mod 1, x E 

for some A G Also write A* G for the adjoint of A. We have for each a; G T'^ that 

k ■ T{x) = k ■ A{x) — k ■ [A(a;)J = A*(k) ■ a; — k ■ [A(a;)J, (4.1) 

where \_y\ G Z'^ of i/ G is the vector with entries \_yj\ for each coordinate j = l,...,d. 
Hence as k G Z'^ the number k ■ [A(a;)J G Z. Now as k G Z*^ and A has integer entries, 
the vector A*(k) G Z*^. Therefore, if we do (4.1) again but now with T{x) in place of a;, 
A*(k) in place of k, we still have that k ■ T‘^{x) is of the form (A*)^(k) ■ x + q for the integer 
q := A*(k) ■ [A(T(a;))J G Z. In general 

k ■ T"(a;) = (A*)"'(k) ■ x + q{x, k, n) 

for some q{x, k, n) G Z for any n G N and x G T'^. This yields in particular that 

exp(27rfk ■ T"(a;)) = exp(27ri(A*)"'(k) ■ a;). 

Let ad = inf{|Aa;| : |a;| = 1} be the smallest singular value of A, which is the same as the 
smallest singular value of A*. Thus we have for any n G N and ?/ G that 


\{ATy\>cTM 


( 4 . 2 ) 
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Then we may estimate 

1 ^ r 

I exp(-27rik ■ {T^{x) - T^(x))) d/x(x) 

m,n=l ^ 

N 

= E / <ixp(-2^!((A-)'”(k) - (A-)”(k)) . x) d^x) 

m,n=l ^ 

^ ^ N m—1 

« E E l«0*)’”(k) - 0-)"(k))|. 

m=2 n=l 

By the mean value theorem, whenever 1 ^ n < m, we can hnd n ^ g ^ m such that 
a™ — — [(^d^ogadlim — n) ^ c{m — n) for c = loga^ > 0 as > 1. Hence by the Fourier 
decay assumption and (4.2) we obtain 


N m—l 


N m—l 


Y, Y l/S((V)”'(k) - (V)”(k))| < ^ |(V)'”(k) - (V)"(k)|- 


m=2 n=l 


m=2 n=l 

N m—l ^ 




N m—l 

EE 


(m-n)“|ck|^ 


where we used |k| ^ 1 and c > 0. Since a > 0 we have that 


/ WT < OO- 
TV 

A=1 

Thus by Lemma 4.1 there is a subsequence Nj —)■ oo, j G N, such that 


E rjv < oo and lim 

^ -7—Vrv 

j = l 


iV, 


i+i 


J^OO 


Nd 


= 1 . 


In particular, the former condition yields that 

« OO 

/ dfi{x) < oo 


i=i 


and so S'at^. (t) —)■ 0 for p almost every a; as j —)■ cx). Now we just check that the latter condition 
actually yields Sn{x) —)■ 0 for p almost every a; as iV —)■ oo, which is what we need. Fix iV G N 
and hnd j G N such that Nj ^ N ^ ^j+i- This yields that 


|TVS'Ar(a;) - NjSNj{x)\ 


N 

exp(-27rfk-T”(a;)) 

n=Nj-\-l 


^N-Nj^ Nj+i - Nj. 
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Hence 

which converges to 0 as j —)■ cx) at almost every x. □ 

5. Further research 

5.1. Gaussian processes. A promising further problem would be to try to extend this result 

for other Gaussian processes. The most probable next step would be to obtain bounds for 
the Fourier dimension of the graph of the fractional Brownian motion = Bf^ for a Hurst 
parameter 0 < < 1. This is a Gaussian process with mean 0 and the increments Bf^ — B^ 

have the variance 

E|Bf-Hf|2 = \t-s\^^. 

Note that in the case H = 1/2 the process Bf^ is just the standard Wiener process Wt- 

Kahane [19] already established an analogue of the bound on the Fourier transform for the 
images of Lebesgue under which means that we could bound the Fourier transform of the 
measure p (Lebesgue measure lifted on the graph G{B^)) for angles 6 which are not too close 
to 0, TT or 271 with |/i(.^)| = A/log |.^|). However, the issue comes from the fact that 

ltd calculus is dehned for ltd drift-diffusion processes for which the diffusion component only 
has the Wiener process. There are analogues of ltd calculus available for Bf^ depending on the 
roughness of B/^, see for example [29] for an overview. 

In the case H ^ 1/2 ltd’s lemma for processes Xf = btdt + atdBf^ has the form 

/(Xt)-/(Xo)= rbtf{X,)dt+ r r atf{Xt)dB/^, 

Jo Jo ^ Jo 

where J Xt dB/^ means a stochastic integral in the sense of pathwise Young integrals, see [38, 39]. 
However, if we recall the proof of Lemma 3.6 when we applied ltd’s lemma, the problem now 
comes that even though bt and at are time-independent and f and / are — 1 or i multiple of 
/, we will have the term in the integral, which prevents us from having an equation 

r f{xt)dt=c r f{xt)dB/^ 

Jo Jo 

for some constant c G C. Burkholder-Davis-Gundy type-inequalities have also been developed 
for fractional Brownian motions, see for example [28]. For small H tools from rough path 
theory might be of help here. Moreover, the recently developed methods by Peres and Sousi 
[30] for the Hausdorff dimension of the graph of fractional Brownian motion with variable drift 
might be relevant. 

5.2. Other measures on the Brownian graph. Theorem 1.3 concerns the push-forward of 
Lebesgue measure on [0,1] onto the graph G{W) but it would be an interesting problem to see 
if one can have similar results for other possibly fractal measures on [0,1]. A possible problem 
could be: 

Problem 5.1. Classify measures r on [0,1] such that for some 0 < s ^ 1 we have 

|f«)l=0(|{|-/2), 1^1 ^oo, 
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such that their lift onto the graph of G(W) under t ha (t, Wt) has power decay with the 
exponent —s/2 for the Fourier transform. 

This is motivated by the fact that in Kahane’s work [19] it is possible to transfer information 
on the Fourier decay (or Frostman properties) of r onto the image measure. Thus for directions 
6 bounded away from 0 and tt we could still bound using Kahane’s work. The main 

problem in generalising our approach to fractal measures r on [0,1] comes from the lack of 
ltd calculus. One would need to hnd a way to link the integrals fj' Zt dT{t) to some form of 
stochastic integrals “jJ' Zt dW/^ involving r. 


5.3. Prevalence. Given that it is often difficult to study the graph of a specihc continuous 
function, one is often interested in the generic case. Indeed, our results on the almost sure 
properties of the Brownian graph are of this nature. When one wants to study generic behaviour 
within the Banach space ((^[0,1], jl-jloo), there are two main approaches: topological (using Baire 
category), or measure theoretic. The measure theoretic approach is often based on prevalence, 
introduced by Hunt, Sauer and Yorke [15], which plays the role of “Lebesgue almost ah” in 
inhnite dimensional Banach spaces in the absence of a useful analogue of Lebesgue measure. 
In our context the dehnition is as follows: a Borel set A C G[0,1] is prevalent if there exists a 
compactly supported Borel probability measure P on Cf), 1] such that 

P(G[0,1]\(A + /)) =0 


for all / G G[0,1]. 

The results of Mauldin and Williams [26] and Orponen with the authors [14] yield that the 
for a Baire generic f & C[0, 1] the graph has Hausdorff dimension 1 and Fourier dimension 0. 
For the prevalence case, Fraser and Hyde [10] proved that the graph of a prevalent / G Cf), 1] 
has Hausdorff dimension 2. This result was later generalised by Bayart and Heurteaux [2], who 
simplihed the approach of Fraser and Hyde by using fractional Brownian motion. The main 
idea is that P is assumed to be the distribution of the fractional Brownian motion t ha in 
(^[0,1] for a given small H, and then one establishes that the graph G{f + B^) of the drifted 
process f{t) + Bf^ has the same dimension as G{B^) almost surely for any hxed / G (^[0,1]. 
By Adler’s results [1], the Hausdorff dimension of G{B^) tends to 2 as iL —)■ 0. 

For Fourier dimension it could be possible to adapt the same approach. Since G(W) has 
the maximal Fourier dimension 1 almost surely, it is likely that that the Fourier dimension 
of a prevalent / G (^[O,!] is 1. The problem reduces to proving that the graph measures 
fVf+w on the drifted process f(t) + Wt almost surely have the same Fourier decay as p, in 
Theorem 1.3 for any hxed / G Gp, 1]. In this case the Fourier transform fif+w has the form 
/i/+w(0 = /o^exp(iYi) dt for the process 

Yt := —27ru{tcos6 + f{t) sin 6* + HA sin6^), 

which is no longer an ltd drift-diffusion process unless f{t) is linear (compare this with (3.1)). 
Thus one would need to hnd ways to develop theory for such processes to adapt the approach 
presented here. 
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